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Dynamics of a polymer in an active and viscoelastic bath 
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We study the dynamics of an ideal polymer chain in a viscoelastic medium and in the presence 
of active forces. The motion of the centre of mass and of individual monomers is calculated. On 
time scales that are comparable to the persistence time of the active forces, monomers can move 
superdiffusively while on larger time scales subdiffusive behaviour occurs. The difference between 
this subdiffusion and that in absence of active forces is quantified. We show that the polymer 
swells in response to active processes and determine how this swelling depends on the viscoelastic 
properties of the environment. Our results are compared to recent experiments on the motion of 
chromosomal loci in bacteria. 


I. INTRODUCTION 

The dynamics of a polymer in a viscous solvent in equi¬ 
librium is well understood mm- Starting from the ex¬ 
actly solvable Rouse model [T], which neglects excluded 
volume and hydrodynamic interactions, it was found that 
the centre of mass of the polymer diffuses with a diffusion 
constant D that is inversely proportional to the number 
of monomers N. The time for the polymer to diffuse 
over its own radius of gyration Rg{^ then intro¬ 

duces a typical time scale, called Rouse timescale, th, 
which scales with iV as ~ ~ iV^. Individual 

monomers subdiffuse for times smaller then r^, but fol¬ 
low the diffusion of the centre of mass after [2] . Going 
beyond the Rouse model, the effects of excluded volume 
and hydrodynamic interactions can be incorporated us¬ 
ing scaling theories [3], simulations and exactly solvable 
models [4]. 

Much less is known about the dynamics of a polymer 
in the crowded [5] and nonequilibrium cellular environ¬ 
ment. The crowdedness introduces viscoelastic behaviour 
[5] with a long term memory while the action of molecu¬ 
lar motors and other ATP-driven, active, processes puts 
the cell out of equilibrium. 

Recently, several experimental and theoretical studies 
have investigated the role of active processes on biopoly¬ 
mer dynamics. On a coarse grained scale, activity can 
be seen as an extra source of randomness in addition to 
that due to thermal motion. This is reflected in a ran¬ 
dom motion of tracer particles (and small biopolymers 
like proteins) that is much enhanced in comparison with 
thermal Brownian motion [THin]- In other experiments, 
the internal dynamics of a polymer was found to be in¬ 
fluenced by the presence of active forces. We mention 
the effect on the bending dynamics of microtubuli m- 
More relevant for the present paper are studies in which it 
was found that the motion of chromosomal loci in simple 
organisms like bacteria and yeast are sensitive to active 
forces [H US]. For example, Weber and collaborators 
[13] found that after addition of chemicals that inhibit 
ATP-synthesis, the diffusion constant of chromosomal 
loci decreased by 49%. Also measurements of chromatin 
dynamics in eukaryotes show evidence for an important 


role played by ATP-dependent processes [13] . 

In the theoretical description of the motion of chromo¬ 
somal loci, it came as somewhat of a surprise that again 
the simple Rouse model turned out to be relevant m- 
Indeed, it was recently argued to give a good description 
of this motion, both between m and during chromoso¬ 
mal segregation m- The reason for this may be the ac¬ 
tion of topoisomerases and related enzymes which cross 
chromosome strands and thus make the bacterial chro¬ 
mosome a phantom chain m- However, other models 
have been introduced to describe the motion of bacterial 
chromosomes, for example, in terms of self-adhesion of 
monomers [18] . None of these models did however ex¬ 
plicitly investigate the effects of active forces. 

In the present paper, we extend recent work on the 
Rouse chain in a viscoelastic medium [nKHHDj by in¬ 
cluding active forces. The advantage of our model is that 
it is exactly solvable and hence like the original Rouse 
model can be used as a benchmark for studies of more 
realistic models that include physical effects like self¬ 
avoidance, bending rigidity and so on. Some of our re¬ 
sults (superdiffusion, polymer swelling) have indeed been 
recently seen in simulations [211 [22] of such models. How¬ 
ever, all the numerical work known to us is for the viscous 
regime. Our work is the first to include both effects of 
viscoelasticity and active forces, two ingredients that are 
necessary for a proper modelling of a cellular environ¬ 
ment. 


II. MODEL 

The Rouse model mm is the starting point in all dis¬ 
cussions of polymer dynamics. It models a long polymer 
chain as a set of N beads (monomers) connected by har¬ 
monic springs. It is important to point out that such 
a description of a real polymer is only appropriate at a 
coarse grained level, i.e. at length scales above the per¬ 
sistence length. Hence, our model will only be relevant 
for long biopolymers (like chromosomes) and also other 
physical properties, like crowdedness, viscoelasticity and 
active forces, will be described at this coarse grained 
level. 

We denote by (f) the position of the n-th monomer 


2 


(n = 0,1,..., TV — 1) at time t. Let us discuss the various 
forces acting on this monomer (for a schematic represen¬ 
tation of our model, see Fig. [^. 



FIG. 1. (Color online) Cartoon of onr model : a large poly¬ 
mer (black spheres connected with springs) in a dense crowded 
cellular environment (orange spheres represent other smaller 
biopolymers like proteins, RNA, ...) that is bronght out of 
equilibrium by active processes (blue spheres with arrows rep¬ 
resent active particles that for a time of order ta move in a 
fixed direction). 

In the Rouse model, the monomers are connected by 
springs with spring constant k. In equilibrium, the av¬ 
erage squared distance 6^ between two monomers then 
follows from the equipartition theorem and equals = 
SksT/k where ks is Boltzmann’s constant and T tem¬ 
perature. A Rouse chain gives a description of a semi- 
flexible polymer on length scales where b corresponds to 
the Kuhn length [23]. We will apply our model to the 
bacterial chromosome, for which the Kuhn length is of 
the order of 100 nm or 300 base pairs (bp) [T7j. Since 
the chromosome of E. Coli has approximately 4.5 x 10® 
bp, the relevant value of N is of the order 10^. 

The cellular environment is a crowded material (am 
[25] in which tracer particles (nano particles, mRNA, ...) 
were found to subdiffusive (for a recent review, see pS]). 
While originally the precise explanation of this behaviour 
was unclear, and various mechanisms where proposed, 
recent experimental evidence shows that viscoelasticity 
and its mathematical description based on a generalized 
fractional Langevin equation best describes the observed 
behaviour [27]. The friction force Fn{t) on a monomer 
in a viscoelastic medium has memory and is commonly 
described in terms of a power law kernel K(t) = (2 — 
a)(l — a)t~°‘ 

Fn{t) = ( 1 ) 

Viscous behaviour is recovered for a = 1 (for which 
K{t) —>• 5{t)) while for a = 0 we get elastic behaviour. 
For 0 < a < 1, we have the viscoelastic situation, inter¬ 
mediate between elastic and viscous response. Estimates 


of a range from a ~ 0.7 for E. Coli m to a « 0.2 in 
the cytoplasm [l^ of eukaryotes. An expression for the 
friction like Eq. Q breaks down on short timescales of 
the order of the molecular collision time scale. In bacte¬ 
rial cells, an upper limit to memory effects is put by the 
cellular lifetime. 

The random thermal force ^T.n{t) acting on the n-th 
monomer is given by a Gaussian random variable, with 
average zero and a correlation that is coupled to the ker¬ 
nel K (t) by the (second) fluctuation-dissipation theorem 

m 

(irAt) ■ A,At')) = AksTKilt - t'\) 5n,m (2) 

With the choice of the power law kernel, Cr.n(i) becomes 
fractional Gaussian noise [29]. 

Much less is known about the precise form of the ac¬ 
tive forces AAt)- f9,ct, the precise characterisation 
of their statistical properties in living cells is a topic of 
current research m l30H33] . Ultimately, the activity is 
due to active ’particles’ (molecular motors, or other ac¬ 
tive proteins) that consume energy to generate motion 
and associated dissipation. In a simple model (see [33] 
for a review), self-propelled particles (SPPs) move with 
constant velocity in a direction e that is subject to rota¬ 
tional Brownian dynamics. This leads to a typical motion 
in which the autocorrelation of e decays exponentially. In 
[22] . the dynamics of a polymer in the presence of such 
SPPs is studied. The polymer and SPPs interact through 
a truncated Lennard-Jones potential. On a more coarse 
grained level, this type of interaction leads to a random 
force on the monomers. The force has an exponential 
correlation with a timescale ta, which can be interpreted 
as the typical time during which the SPPs move in a 
straight line. 

The modelling of active processes through this type of 
random forcing is quite common in the literature and was 
for example used in m in simulations of active semiflex- 
ible polymers. We will follow this approach since our 
model, as already stressed above, is defined on a coarse 
grained scale. In conclusion then, we assume that the 
active force on the n-th monomer, is a Gaus¬ 

sian random variable with average zero and a correlation 
given by 

(AAt) -AAt')) =iCe^^{-\t-t'\/TA)5n,ni (3) 

Here C characterises the strength of the active noise. We 
do not include possible spatial correlations in order to 
keep the model soluble. Moreover, little experimental 
insight on such correlations is available. 

Putting everything together and neglecting inertial 
terms, the equation of motion of the n-th monomer is 
the overdamped generalized Langevin equation 

FAt) - k {2Rn{t) - Rn+i{t) - Rn-i{t)^ 

F AAt) F A,nit)F[{t) = 0 (4) 

where F[(t) is the Heaviside function. It is important to 
point out that the active forces, since they are not related 
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to the friction kernel, put the system out of equilibrium. 
Starting in equilibrium at t = 0, the solution to equation 
@ therefore gives the response of the polymer to, for 
example, the addition of ATP at t = 0. After a long 
time, the polymer will evolve to a new, nonequilibrium 
steady state. In the next section we will calculate both 
the transient and steady state behaviour of the polymer 
after activation of the active forces. 


III. RESULTS 

The techniques to solve the set of equations 0 with 
appropriate boundary conditions are standard [2j. The 
details are given in the supplemental material to this pa¬ 
per [ 33 . Here, we discuss only the results. 

In the dynamics of the polymer two relevant time 
scales occur. The first is ta, the second is the 
Rouse time t/j, which in the viscoelastic medium equals 
— a)jb'^N'^ [THl EH] . Notice that for 
the bacterial chromosome, using the values of b, N and a 
quoted above, tji can become quite large, indeed larger 
than the duration of the cell cycle. The persistence time 
of active processes should be much smaller then the cel¬ 
lular lifetime, so that when applying our model to the 
bacterial chromosome, the two timescales are separated, 
Ta ^ Tn. Experiments on chromosomal loci are done on 
timescales of O.I seconds to minutes, i.e. in the regime 
~ t < Tr. 

We first discuss the motion of the centre of mass 
Rem (t) ■ It was found earlier that in a viscoelastic medium 
but in absence of active forces m, the centre of mass 
Remit) performs a subdiffusion 

aUt) = ( - RUO))") = (5) 

where Ga = r(Q;)r(3 — a). To this subdiffusion an extra 
term is added when the active forces are turned on. It 
equals 



where r(Q!; y, x) is a difference of two incomplete gamma- 
functions. Firstly, observe that as in the standard Rouse 
chain, the (generalized) diffusion constant remains in¬ 
versely proportional to N. More interestingly, it can be 
shown that for t <C (|^ goes as i.e. is superdif- 
fusive if a > 1/2. On the other hand, for t ^ ta, ®> 
evolves as i.e. slower than ([^, so that asymptot¬ 
ically in time, ^ t“. The resulting behaviour is 

shown in Fig. for a = 0.7, C = 100 and various val¬ 
ues of ta [36]. As the figure shows (boxed region), if 


Cta is large compared to over several orders of 

magnitude in time, the centre of mass subdiffuses with 
an exponent 2q; — 1, whereas without active forces, the 
exponent would be a. If this time regime corresponds 
to the experimental one, not taking into account active 



FIG. 2. (Color online) Log-log plot of squared distance 
travelled by the centre of mass as a function of t/Tjj for a 
Rouse chain in a viscoelastic medium (a = 0.7, ksT = 7 = 1 ) 
in the presence of active forces with C = 100 and taItr = 
10“^, 10“®, 10“"^ (full lines, top to bottom) compared to that 
without active forces (dashed line). The inset shows the same 
for the viscous case, a = 1. Results are for N = 256. 

forces could lead to a wrong estimate of the exponent a, 
and hence to a wrong characterisation of the rheological 
properties of the environment. Results like that of Fig. 

are under the assumption that our model is realistic 
for all t. When comparing with experimental data, one 
has to take into account that, as already discussed in 
the previous section, the model will break down on small 
and large time scales. Therefore in the cellular context, 
it is possible that only the superdiffusive behaviour is 
observed. 

For the viscous case, 2a — I = a, so that after an 
initial regime of ballistic motion cr/^(t) ^ the polymer 
performs ordinary diffusion but with a diffusion constant 
that is enhanced by a factor I -|- (see inset of 

Fig. [|). 

A second global quantity is the end-to-end vector 
P{t) = Roit) — RN-iit). Its averaged squared length, 
R^{N,t) = {Pit) ■ Pit)), measures the size squared of 
the polymer and equals in equilibrium, both in viscous 
and viscoelastic media, b^N. In response to active forces, 
R^ {N, t) gets an additional term which equals 
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FIG. 3. (Color online) Squared end-to-end distance of the 
polymer as a function of time. In the main figure (inset), the 
equilibrium length is (not) subtracted. The straight line has 
slope 1.4. {ta/tr = 10"'^, C = 100, a = 0.7, N = 256) 


Here Ea^p{z) is the generalized Mittag-Leffler function 
m, and Tp = Since 0 is positive, we conclude 

that active forces swell the pmymer. It can be easily 
shown that initially this swelling grows projrortional to 
after which R^{N,t) saturates. In Fig. ^ we show 
the results of a numerical evaluation of ([^ for a polymer 
with N = 256 in a medium with a = 0.7. 

The value at which R^(N, t) saturates, i.e. the squared 
size of the polymer in the new, nonequilibrium steady 
state, has an interesting TV-dependence. For the viscous 
case, a = 1, the integrals in 0 can be easily calculated. 
In this way, it is found that the swelling of the polymer 
is proportional to N so that 


RUN) = ^hm i?2(iV,t) = (^1 + b^N (8) 

In fact this result is exactly what can be expected from 
a recent extension of the equipartition theorem to har¬ 
monic oscillators in viscous, active media [38) . In that 
reference, the average potential energy of a harmonic os¬ 
cillator in a viscous, active bath is calculated and com¬ 
pared with experiments. Using the results of that paper, 
and the fact that the Rouse chain consists of N indepen¬ 
dent harmonic oscillators one can also derive (|^. 

In the viscoelastic case, the situation is more compli¬ 
cated. It is possible to determine the leading behaviour 
of the integrals in Q for TV 1 (and t —)■ oo). In this 


l^/a-l 







( 7 ) 


way it is found that for 2/3 < a < 1 


RUN)=b^N + 


48(7rA(47r^fc)^/' 
(77(3 — a))i/° 


a-2 


(9) 


where 


poo 

f{a) = C(4 - 2/a, 1/2) / dx x‘^°‘-‘^EU 

JO 

and C,{x,y) is the Hurwitz zeta function. In Fig. 13 
we show the results of a numerical evaluation of (7) 
for a = 0.8 together with the asymptotic behaviour (9) 
(dashed line). For a < 2/3, numerical integration of (7) 
indicates that the swelling of the polymer approaches a 
constant as TV increases (Fig[^. From our calculations, 
we see that in an active environment, large polymers 
are orders of magnitude more compact in a viscoelas¬ 
tic medium with small a. This observation could be of 
relevance for storing a large chromosome in a small cell. 



N 


FIG. 4. (Golor online) Difference between squared end- 
to-end distance of the polymer in the nonequilibrium steady 
state and in equilibrium as a function of TV for a = 1,0.8 
and 0.4 (top to bottom). The symbols are the result of a 
numerical evaluation of (ItI for T —^ 00 . The full line is a plot 
of 1^, the dashed line of (|^. 

Finally, we turn to the motion of the individual 
monomers. From an experimental point of view, this 
quantity is the most interesting one, since it can be de¬ 
termined using fluorescence techniques na na [391 HO]. 
We present results for the middle monomer, but the be¬ 
haviour of other monomers is qualitatively the same (see 
supplemental material for the general case |35j). 
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FIG. 5. (Color online) Log-log plot of o'm(t) as a function of 
time in presence of active forces (full lines) for a = 0.7, ksT = 
7 = 1,(7 = 10'*, iV = 256 and ta/tr = 10■^ 10"^ lO^* (top 
to bottom). The dashed line is the same quantity in absence 
of active forces. The inset shows that the results for a smaller 
time window in absence and presence of active forces {C = 
1,ta/tr = 10“^) run parallel. 


In absence of active forces, one has m 


^mit) — ( {RN/2{t) - ^Ar/2(0)) 


46^ iV ^ J_ 

TT^ 

p—2,even 


l-E, 


a,l 


-fNaGa 
t 


( 10 ) 


For t < tr, the monomer is found to subdiffusive with 
an exponent a/2 whereas on time scales larger than 
the monomer follows the motion of the centre of mass, 
hence subdiffuses with an exponent a. Since for the cy¬ 
toplasm of E. Col/ a equals approximately 0.7, the ob¬ 
served exponent of the subdiffusion of the chromosomal 
loci, 0.39 ± 0.04 [131 m], is consistent with a/2 within 
the experimental error. 

In the presence of active forces two terms have to be 
added to ([To| . The first one is the same as for the centre 
of mass the second equals 1/4 of Q (but now with 
the sum over the even modes). 

Analysing the resulting behaviour, we find that as was 
the case for the centre of mass, the short time response of 
the monomers to active forces is superdiffusive (at least 
when a > 1/2). This behaviour, where ~ 

holds for t < TA- Very recently superdiffusive motion of 
chromosomal loci has indeed been observed [53] ■ It is 
not clear whether these rapid chromosomal movements 
are due to active processes or to stress relaxation. Our 
results however quantify better the response to active 
forces and could therefore help in discriminating the real 
origin of the observed motion. 


For t ^ tr the monomers follow the subdiffusion 
of the centre of mass. The behaviour in the intermedi¬ 
ate time regime, ta ^ t ^ tr, is more complicated and 
contains terms proportional to (coming form ((lol)), 
^ 2 a-i centre of mass), t3a/2 ^3a/2-i_ 

simple power law behaviour emerges. A plot of the full 
expression for cr^(t) (Fig. for C large shows that the 
initial superdiffusion is followed by a behaviour, af¬ 
ter which there is a large time window in which the 
term dominates. The crossover time to the ^“-regime is 
very large, and may not be observable for chromosomal 
loci. For certain values of C and ta, there is a time regime 
(see inset of Fig. where, on a log-log scale, the graphs 
of cr^ (f) in presence and absence of active forces run par¬ 
allel. This scenario resembles, at least qualitatively, the 
experimental one, where after inhibition of ATP synthe¬ 
sis with sodium azide and 2-deoxyglucose, the exponent 
of the loci’s subdiffusion hardly changed but the diffusion 
coefficient decreased m- 


IV. CONCLUSIONS 


In summary, we have studied the behaviour of a long 
ideal polymer chain in a viscoelastic and active bath. We 
have formulated and solved a nonequilibrium version of 
the Rouse model, a model that in equilibrium forms the 
starting point of the theory of polymer dynamics. The 
results obtained are interesting for nonequilibrium sta¬ 
tistical mechanics (power law response to active forces, 
equipartition theorem out of equilibrium) and polymer 
physics (swelling of polymer in an active bath). More¬ 
over, our results show a qualitative similarity with exper¬ 
imental results on bacterial chromosomes. We mention 
the regime with superdiffusive motion of monomers and 
the fact that in certain time frames the exponent of the 
subdiffusive motion is almost independent of the pres¬ 
ence/absence of active forces. 

While the Rouse model has been claimed to describe 
well several properties of bacterial chromosomes [muz], 
it is clear that such a description cannot be valid for all 
time and length scales and that effects of self-avoidance, 
bending rigidity, topology and so on have to be taken 
into account. This is even more true for chromatin for 
which modified Rouse models have been introduced that 
take into account bending rigidity [H] or long range 
interactions m- The motion of chromosomal loci was 
also investigated on the basis of scaling arguments and 
computer simulations [ 32 ] in the fractal globule model, 
which shares topological properties with chromatin [43] . 
To investigate these issues further, we have recently 
developed an algorithm that allows to simulate particles 
and polymers in a viscoelastic medium subject to 
correlated noises coming both from thermal forces and 
active forces. With that algorithm it becomes possible to 
include effects of self-avoidance and bending rigidity and 
to see how the behaviour found here for an ideal chain is 
modified. The results will be published elsewhere [33] ■ 
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